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1.  Introduction 

Consider  the  following  elliptic  boundary  problem 

|  ~  ~ rt*’*)'  (**»)€  Dm  {(*,*)  |°  <*,*,<  1};  ^  ^ 

1  «»(*,  y)  -  fix,  if),  (*,  y )  €  SD-, 

where  p,-(*,y)  >  po  >  0,  «  =  1,2  and  po  is  a  constant.  Approximating  (1.1)  by  a  five-point 
finite-difference  scheme  on  the  grid  system  Du  m  {(jA,  kh)\  h  «  l/(n  +  1);  j,  k  »  1, 2, . . . , »} , 
we  obtain  a  set  of  hear  algebric  equations 


Au  -  Q,  (1.2) 

where  u  and  Q  are  n*-cohimn  vectors  whose  components  are  grid  functions  Ujj,  and  re¬ 
spectively,  ordered  in  a  left-to- right  and  down-to-up  fashion,  and  the  matrix  A  is  defined  by  the 
following  * 


(An)/,*  ■  S/,*n/,*-i  +  DjjkUj-xj,  +  £/,*«/,*  +  #/,*«/+M  +  (1-3) 


whose  coefficients  satisfy  the  conditions 


(l)fi/*,  %,  F/,*  and  G/,*  are  all  less  than  a  negative  constant  if 
subscripts  of  their  associated  grid  functions  are  in  [l,n]; 


(2)£y>  +  BiJt  +  D,,.  +  fy.,  +  Git  { 1 }■ 


BiJ>  •  %*  ‘  F>\*  •  G>,*  *  0 
otherwise; 


(1.4) 


and 

0/+i,*  ■  fy,*»  Bi,k+i  ■ G/,*-  (1.5) 

It  is  well  known  that  Ana  symmetric  positive-definite  matrix  with  five  nonzero  diagonals. 

In  [4],  H.  L.  Stone  proposed  an  iterative  method.  Strongly  Implicit  Procedure  (SIP),  for  solving 
(1.2).  First,  he  approximately  factored  A  as  a  product  of  a  lower  triangular  matrix  I  and  an  upper 
triangular  matrix  V  which  have  tbe  same  sparsity  as  the  lower  and  the  upper  triangular  parts  of 
A,  respectively,  i.e. 


(£«)/,* 

(u")jjk  -»**  +  «yA«y+M  +  /y>*»y,*+i. 


(1.6) 


so  that 

((Ilf  My,*  ■  *y,*«y>-i  «/+,,*-i  +  «/>«/-«.* 

+  (*y.* + *y^/y>-i  +  cy,*«y-i,*)«y,* 

+  d/^e/4U/+i>  +  +  dyA/yA*y>H- 

Let  the  matrix  B  be  defined  as  fbQows: 

(£«)/,*  *  *y^y>-i  {«/+»>-»  +  «»yA  -  «»y+iA  -  } 

1  f  1  ,  \S},kZn 

+  e/.*//-i.*{«/-i>fi  +  «t/A  -  ®«y-w  ~  <"»y,*+i)- 


‘Tlngtnt  «Ut  sepw,  a*  last  as  a  an Ufa  h  Mud  is  tils  say  at  ante  a  eaavttoa  that  if  tha  nbacripta  of  a  grid 
ftnctiaa  aaa  aat  is  fl.s|  this  tha  tanaa  whteh  lariats  aach  pM  lhaatlaa  4a  sat  ariat,  or  coagrimta  art  ana. 
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The  nonzero  elements  of  L  and  V  are  defined  so  that  A  +  B  —  LU.  The  algorithm  for  computing 
these  elements  is 


*j\*  -  Bj,h  -  a*j,* 

«/,*  »  Djjk  -  acj'kfj-ij,  f 

*j,k  +  lj.kfi.h- 1  +  cj.kti-i,k  ■  Ej.h  +  a*i,*«i,*-i  +  aej.kfj-i.k,  l<j,k<n,ae  (0,  lj  (1.8) 

d, .*/,.*  »  <?/,*  -  *Cjjkfj-lJ>, 

The  SIP  iterative  scheme  is  given  by 

(A  +  -  «<**)  -  r(Q  -  y|«W)t  ««0,1,  •.  (1.9) 

The  practical  numerical  computations  given  in  (4]  suggested  that  SIP  is  convergent  and  more  effi¬ 
cient  than  the  standard  iterative  methods.  After  Stone  the  theoretical  analyses  of  the  approximate 
factorization  iterative  methods  given  by  T.  Dupont,  R.  Kendall,  B.  Rachford(2],  A.  Bracha- Barak 
and  P.  Saylorfl]  etc.  were  important  contributions  to  the  later  development  of  the  more  efficient 
preconditioned  conjugate  gradient.  However,  these  studies  did  not  prove  that  SIP  is  convergent 
and  SIP  was  layed  aside.  The  purpose  of  this  paper  is  to  analyze  the  convergence  of  SIP  by  the 
methods  of  those  mentioned  above.  In  Section  2  we  discuss  a  property  of  A  +  B,  where  A  is  not 
confined  to  be  a  symmetric  matrix.  In  Section  S  the  convergence  of  SIP  is  proved. 

1.  The  property  of  A  +  B 

As  in  Lemma  1  in  (1]  we  hare  the  following  theorem. 

Theorem  2-1.  If  the  matrix  A  satisfies  (1.4)  and  0$a  ^  1  then  the  quantities  fay,*,  ey,*,  dy,*,  ey,*, 
and  fj'k  defined  fay  (1.8)  satisfy 

(1)  *y.a  S  Bjjk  and  ey,*  £  £>/,*, 

(2)  dy,*  >0, 

(3)  -  1  <  *yf*  $  FjjkHi*  -  1  <  fi.h'Si  Gj.kMjJt, 

(4)  1  *f  ey,*  +  fjjt  >  0, 

and  they  are  all  cootmaom  faactictm  of  a  €(0,1  j. 

Proof.  The  proof  is  by  induction.  Obviously,  the  theorem  holds  for  /  »  fa  •  1.  Now  assume  the 
conclusions  of  the  theorem  hold  for  all  points  preceding  (/,  fa).  First,  we  have 

*i4  ■  By,*/(1  +  £  By,*, 

*  By^/(1  +  afj-\,h)  S  By,*, 

and 

*j+  mEjJ>~  *j,kfjji-i  -  CjAtj-i*  +  o*y,*ey,»-i  +  **jjhfj-\jk 

“  By>  -  hjjkfjjt-x  -  +  (By,*  -  fay,*)  +  (By,*  -  ey,*)  (2.2) 

■  By,*  +  Bjj,  +  By,*  -  *y>(l  +  //,*-»)  -  ey,*(l  +  «y-t,*)  >  0. 
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Since 


*j*  +  dy.*«y.*  -  By,*  +  By,*  +  Dj*  -  *,\*(1  +  /y,*-t) 

-  */,*(!  +  «y-»,*)  +  By,*  “  **y,*e/.*-i 
■  By,*  +  By,*  +  By,*  +  Fy,* 

-  *i.*(l  +  Qe/,*-»  +  //,*-») 

“  ei,*(l  +  ey-i,*)  >  0 

then 

1  +  eyjk  >  0, 
i.e. 

-1  <  tjj,  »  (Fy,*  -  Q*y,iey,*-i)/rfy,*  5  Fy,*/dy,*. 

Similarly, 

—l  <  /y,*  — 

In  the  above  expressions,  both  numerators  and  denominators  are  continuous  functions  of  a  € 
[0,1]  and  the  denominators  are  nonzero.  So  ty,*,  cy,*,  rfy,*,  ey,*  and  /y,*  are  all  eontimioas  functions 
of  a  €  [0,1]. 

Finally,  since 


rfy.*(l  +  «y.*  +  ft*)  *  ^y.*  +  rfy,*«y>  +  dy,*/y,* 

■  JEy,*  +  By,*  +  By,*  +  Fy,*  +  Gy,* 

“  *y,*(l  ♦  Qei>-i  +  /y,*-i) 

“  fy,*(l  +  *i-\*  +  tt/y-t,*)  >  0 


then 

The  theorem  is  proved. 


1  +  «y,*  +  /y,*  >  0. 
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Theorem  M.  if  the  matrix  A  satisfies  (1.4)  and 

By,*  +  By,*+1  +  By+,,*  +  Fy_if*  +  Gy,*_,  ^  0, 

( i.e.,  Ar  is  diagonally  dominant;)  1  <;’,i  Sn  (2.3) 

By,**.,  >  Gy,*,  By+i,*  >  By,* 

and  the  matrix  B  it  defined  by  (1.7)  and  (1.8),  then  there  exists  a*  €  (0,1)  mch  that  fora  €  [0,a*j 
the  matrix  A  +  B  is  positive  real  mid  satisfies 


((A  +  B)*,*)  >  (1  -  fi)(Am,«)  (2.4) 

fcr  some  fi  satistyinf  0  £  <  1  and  aay  real  vector  a  ft  0. 

Note  that  (2.3)  is  weaker  assompthm  than  (1.8)  but  (1.5)  impBee  (2.3). 
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Proof.  For  any  real  vector  u  ft  0  we  hare 


(y4w, «)  »  +  5Z5Z  DiW-i.W  * 

/■l  t>>  />>  Imi 

It— 1  It  t  ■— 1 

+  EE  +  EH  Gitkoj'k+lmi<i 

jml  M  /at  M 

+EEM,* 

jml  km  | 

*  “5  E E  si.*{(“y>-»  -  «M)a  -•£*-  «Lt-» } 

jml  km* 

-  5  E  E  -  «/,*  -  «>-i,*} 

/aj  Sal 

“  5  E  E - “/.*)* “ «J,*  - «*+m} 

jml  kml 

-  |EE°w{^**w  -  */>)*  -4*-  «i.*+i} 


jml  kml 

;«1  M 


Making  some  transformations  of  subscripts  in  the  above  expression  we  get 


Mw, «)  -  -  5  2  +  Cy,*)(«y*M  -  «>,*)* 

/at  Sat 
«-l  * 


“  5  ]C  D***m  ♦  -  «y,*)* 

jmkrni 

+  5  E  EWw + Bi* + + + 

y- 1^» 

I  a  • 

+  2  E  D1**  +  fiy>+» +  ^V+m  +  ^y-14  +  Gy,*-  i)«l*  >  °- 


jml  Sal 


(2.5) 
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On  the  other  hand,  for  any  real  vector  ti  j*  0  from  (1.7)  we  have 


a-i  a 

(Bu,u)  m  53  E***^-*”1 {“/+»>-»  +  avS,l>  ”  «*»/+!>  ~  1 }“/,* 


j»  1  M 


*  r  ■» 

+  E  E  C/.*//-m{u/-1,*+1  +  ««/,*  -  -  ««/,*+!  }“y.A 

/at  tal 

l  *_l  *  r 

-  -  5  EE  “  «/.*)* 

yal  teat 

-  q («>,*-,  -  «>,*)*  -  «£*  -  +  ®«*+u  + 

l  *  *”1  r 

-  2  5Z  E  [(«/.*  -  )*  -  «(«/-!,*  -  “/>)* 

y«*  **t 

-  «(«y,*+i  “  “/,*)*  ~  “y.4  -  “y-^te+i  +  «$-M  +  <*«y>n]  • 


Varying  the  subscripts  and  using  the  convention  *y,i  and  cij,  »  0  we  get 

I  n-I 

(Bo, a)  -  -  -  53 E(*y>M«i.A +«y+i,*/y,*)(«*y,A+i  -  «y+i,*)* 
y«i  m 

+  2  E E^****-1 + ey+i.*/y^)(“y+i4  -  «y,*)z 
4  y»i  M 
a— I 

+  1 53  *y,a«y^-i(«y+i,»  -  *,>)* 

y-» 

a  t*l  *”* 

+  »  E5 Z(*/.*+i«y.»  +  cy,*/y-a)(«y>M  -  «y,*)z 

y*i  tei 

Cl  t{ 

+  n"  E  e*i*/»“l.*(*%*+l  _  **•»*)* 
tel 

+  5  E  E  *y>ey.*-i« y^  +  5  E  E  *y-i>+i*y-M«y> 


y-ites 

•— 1  a— 1 


2 

y«*  tei 


-  ? EE *y.tei«y4«J> “  2  EE *y-M«y-M-i«y^ 

y-i  *-i  y-stes 

+  5  E  E  ey4/y-i.*«y,A  + « E  E  ey+»*-t/y.*-»«y,A 

y-*  *-i  y-i*-* 

a— la— 1  •  » 

-  2  EE «y+u/y^^ “  o  EE ey^-i/y-»>-i*?.A- 

yal  tel  yatftes 
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Dividing  each  of  the  laet  right 
for  example 


latione  of  the  above  expression  into  four  parte,  respectively, 


1^»  »  |  a-la-I  *-l 

2  */,Ae/, *-!«/,*  ■  2  {£ H i«£*  +  £ 


2 

/si  te3 


jmtkm t 


»-*l 


and  defining  Sj,  Sj,  Sj  by 


I  «-i  ■-! 

(5,u,u)  z  -  -  +  e/H>//>)(«/,H»  ~  «/+ 1.*)* 

*  /si  tel 
^  h"I  •— I 

+  2  5Z5Z^**m->  +  c/+u/m)(«/+M  -  «/,i)J 

*  /•»  tei 

+  2  51  2](*/,Hl«/,A  +  C/.*/i-U )(«*/,*+!  -  «/,*)* 

/si  tel 

H  «-l 

5J  (•/+»,»  ~  «/,■)*  +  «■  c%*/»-ijk(vfi  ~  *»,*)2 

/■»  tel 

+  +  Ct,»-1  /l,»-l )*?,»  +  |(e«,l/«-M  +  V-W«a-l,l)<i, 


1 

(5a«,u)  z-  53  51  [(*/.*e/^-i  +  */-i,*-He/-i>)  -  o(*/-M«/-i,*-i  +  */, Hi </,*)]«/,* 

/tetes  J 

I  *-l  ■— 1 . 

+  2  51  £[(*/•*/>'-».*  +  ei+i,*-j/»>»)  ~  ®(«/+m/i\*  +  e/,*-t//-i,*-i)l  «/,* 

/■J  tet 

i  r~i  r 

+  2  5Z  [(*i-»4e/-M  -  «*/4e/.i)  +  (c/,i//-u  -  <m/+m//,*)]u/,i 

/s»  J 

1  1  r 

+  2  a  “  ®*/-i,«e/-i^»-i)  +  (e/+i,»-i//,»-i  ~  «/,■»-! //-i,i»-i  )]•*/,» 

M 

1 r 

+  2  5Z  ”  ®*1>+i«m)  +  )]«?,* 

1  Sir 

+  2  »>♦•»•*-»>  *  ®K-i4«»-i>-i)  +  («m/*-i,*  ”  oc»4-i/--i,*-i  )]<. 


we  have 


(Bn,  u)  -  (5j«,  h)  +  (St u,  u)  +  (S»u,  u) .  (2.10) 

Obviously,  (Sju, a)  >  0.  From  Theorem  2.1  we  know  that  for  the  grid  point  (j,k),  2  <  j,k  <  n-1, 

,  ,  t  ^  Bj,kFj,k-i  .  Bj-t  k+i  Fy_i  * 

*/>«/.*- 1  +  + - J/7^ - >  °- 

This  expression  is  a  continuous  function  of  a  €  (0, 1],  so  it  has  a  positive  minimum.  On  the  other 
hand,  *y_i, *«,_!,*_!  +  1/4+ie/,*  is  also  a  continuous  function  of  a  €  [0, 1].  Therefore,  there  exists 
some  ay,*  €  (0, 1)  such  that  for  a  €  [0,ay,*]  we  have 

(*/4e/,*-i  +  iy-i,*+i«y-i,*)  -  a  (iy-i,*«y-M-i  +  */,  *+!*/,*)  £  0. 

Taking  ai  *  min«<y,*£n-i  {ay,*},  if  a  6  [0,at],  the  first  summation  of  (S*ti,ti)  is  not  less  than 
zero.  Dealing  with  the  other  summations  of  (5sti,  u)  in  the  same  way,  we  can  get  02  €  (0, 1)  such 
that  when  0  <  a  <  <*2,  then 

(Ss«,«)>0.  (2.11) 

In  order  to  discuss  (Siti,  u),  we  cite  Lemma  3  in  [2]  here:  let  c  and  /  be  positive  and  let  a,  I, 
e  be  complex,  then 

^la-kpS/la-ef  +  elk-el*. 

So,  for  each  grid  point  (j,k),  1  <  j,k  £  n  -  1,  we  have 

1  *y,*+t  +  <fy,*ey,*  dy,*«y,*  •  *y,*+*  x2 

2  5^ 

*  -  5  («f*«  -  «/,*)*  +  d/,*e/.*(«*/+M  -  »/,*)*]  (2.12) 

“  ~\kf  k+1  j.dJ'ie*’k  "  «*/,*+!«/, *)(«/,*+!  -  “/,*)* 

+  (*>.*  -  a*/4e/,*-i)(“/+i,*  -  ay,*)*]  • 

Bat,  by  assumption  (2.3)  we  have 

+  *>**}*  +  *#>+>  m*i*  +  +  ^y,*+t  -  ®*y^+i«y,* 

*  d/4  +  d/4«/,t  +  C/4  -  «*y> f iey,* 

*  d/.*!1  +  •/,*  +  //,*)  +  “*/, *//-»,*  -  akjMltjih. 


Because  1  +  e/4  +  //,*  is  a  continuous  function  of  a  €  (0,  Ij  and  haa  a  positive  minimum,  as  in  the 
above  djaetmaion,  there  exists  ay.*€  (0, 1)  anch  that  when  0  £  a  £  a*  *  we  have 

d/4  +  +  */>n  >  0 


-  0/.A  <  1. 


(2.13) 
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Set 


*» 


then,  for  0  <  a  £  a$  we  get 


max 


M 


<  i 


(2.14) 


and 


.  »->»-! 

”5  IZ IZ */.*+»«;•> («;+>.*  -  «/>h)* 

1  jmt I  tel 

.  a— 1  a— I 

t  y.  r\B^  -  oi;,*+i«y>)(«»+i  -  «*/,*)* 

/at  M 


.  a— I »— 1 

+  1 2 “  «*/>«/.*-i  )(«/+!.*  “  «/.*)*• 

Z  /at  tel 


(2.15) 


In  the  same  way  we  can  prove  that  there  exists  <>4  €  (0, 1)  such  that  when  a  €  (0,  «<)  we  have 


o<A<i 


j  a-l  a— 1 

-  o  mz  c/+i,*//.*(u/+‘.*  -  m/>+i)j 

1  /at  tel 

-  5^*  -  oc/+i,*/>\*)(“/+t.*  ~  •/.*)*  (J-16) 

/at  tel 

+  \fo  £  -  «/,*)*• 

/«I  tel 


Let  ^  a  max(0ulh)-  I*  «•  obvious  that  there  exists  as  €  (0,  l)  sneh  that  for  a  €  [0,a*l  we  have 


+  Ft,  1  +  0M)  -  5®(*»4«w  +  e*4/u )  *  0. 

^4(£<M>  +  4m  +  Am»)  "  j®(V-M*a-M-l  +  'a,a-l/e-I,a-l)  ^  0* 


(2.17) 


t 


Finally,  take  a *  ■  min  {a*,  as,  04,05).  Then,  if  0  £  a  ^  a*,  we  have 


1  *~i  *-*t 

(Au,u)  +  (B«,n)  >  -  t(I  -  $)  £  +  G/\*)(uy.*+i  "  «*;,*)* 

i»i  tei 

“  5O  **  P)  £  +  f/a)(*/+M  “  *»/>)* 

y»i  a«i 

1  "-l 

~  r  £(Ba>H  +  Ga,i)(u>»,*+1  ~  «»,*)* 

tel 

1  *_l 

“  5  +  ^/.a)(“/+i,a  ** 

#-» 

j  a— 1  a— I 

+  «o(i  -  If)  £)  +  ej+i.*/i,*)(«y+i,*  -  «/.*)* 

y*i  tei 

.  a— 1  a— l 

+  -0(1  -  0)  £  +  e/>//-i.*)(«/.tei  “  «•/>)* 

y«i  tet 

+  5(1  -^)Z!  +  */.*  +  Di'k  +  + 

y«i  tei 

+  =  Z3  +  %*+* +  9»>*a  +  F>'-*-‘ + Gi.*-»)“yjk 

>-»  *-> 

+  (**,«)  +  (**,«) 

i  (1  -  $){Au,v)  >  0. 


(2.18) 
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S.  Cmnt(m  eeadlttoas  of  SIP 

Now  we  retan  to  the  case  that  the  matrix  A  satisfies  the  conditions  (1.4)  and  (1.5),  so  it  is  a 
symmetric  positive  definite  matrix.  First  we  qpote  Lemma  2.1  in  (3]  as  one  of  our  theorems. 

Theorem  Id.  The  iteratm  scheme  (1.9)  i  ccurer*ent  if  and  only  if 

Be  + *))}>£,  (3.1) 

where  Xj(W)  dbnoOas  m  dgsoidbs  of  the  matrix  W,  and  Bsp  denotes  the  real  part  of  the  complex 
#*• 


Theorem  M.  if  the  matrix  A  setisfiss  (1.4),  (1.5)  end  the  matrix  B  is  defined  by  (1.7)  and  (1.8), 
then  there  exist  •**  <  (0,1)  aad  5**  €  (0, 1)  such  that,  for  0  £  a  £  a**, 

Be  {A/(A“*(A  +  B))}2(1-  r  )•  (82) 


Proof.  As  in  the  proof  of  Theorem  2.2,  we  can  get  a"  €  (0, 1)  and  £**  G  [0, 1)  .  If  0  <  a  <  a**, 


((A  +  B)a,  «)>(!-  fi*)(Au,  u). 


(3.3) 


where  u  /  0  b  any  real  vector.  Let  Ay  be  an  eigenvalue  of  A-1  (A  +  B)  and  *  +  if  ?  0  be  the 
associated  eigenvector,  i.e. 

A“‘(A  +  B)(x  +  if)  »  Ay(x  +  if). 

So 

((A  +  B)(*  +  if),  *  +  if )  ®  Ay  (A{x  +  if),  *  +  if) . 

Because  A  b  symmetric  positive  definite, 

((A+ B)x,x)  +  ((A  +  B)y,y)  *  ReXj{(Ax,x)  +  (Ay,y)} . 


From  (3.3),  we  get 


I 


Based  on  Theorems  3.1  and  3.2  our  main  result  b  obtained. 

Theorem  U.  If  the  matrix  A  satisfies  (1.4),  (1.5)  and  the  matrix  B  is  defined  by  (1.7)  and  (1.8), 
then  exist  a**  €  (0, 1)  and  $**  €  (0, 1).  When  0  <  a  <  a**  and 

r  <  2(1  -  /T),  (3.4) 


the  iterative  scheme  (1.9)  converges. 
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